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We show theoretically how quantum interference between linearly coupled modes with weak local 
nonlinearity allows the generation of continuous variable entanglement. By solving the quantum 
master equation for the density matrix, we show how the entanglement survives realistic levels of 
pure dephasing. The generation mechanism forms a new paradigm for entanglement generation in 
arrays of coupled quantum modes. 
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Entanglement is a key concept in quantum physics and 
is a crucial resource for quantum information science, 
particularly within recent schemes based on initial multi- 
partite entangled states [J, [^j • The generation of entan- 
gled states of two or more quantum modes typically relies 
on parametric down-conversion in nonlinear crystals [3- 
0] or optical frequency combs ■ Schemes working at the 
microscopic scale |8l4ll| - suitable for integrated devices 
- are instead always based on the cascaded biexciton- 
exciton radiative decay in semiconductor nanostructures, 
and thus restricted to bipartite entanglement. 

In systems of weakly nonlinear coupled quantum 
modes, the interaction energy associated with two quanta 
is smaller than the broadening introduced by the finite 
lifetime of the mode. The opposite situation has re- 
cently been the object of theoretical investigation, be- 
cause of the possibility of engineering strongly correlated 
many-particle states, and numerous applications ranging 
from the photon blockade effect |ld 4l4l| to the perspec- 
tive of a quantum simulator [15l - ll8| . The requirements 
for a practical realisation of such a strong nonlinearity 
within a solid state technology are however very strin- 
gent, and perhaps the only clear-cut observation of the 
photon blockade has been reported in a state-of-the-art 
atomic system [l^l . We have recently suggested that pho- 
tons with strongly sub-poissonian statistics can be emit- 
ted by a set of coupled modes in the weakly nonlinear 
regime [l^, thanks to the interplay of the weak non- 
linearity and quantum interference [20[. We argue that 
the same mechanism can be more generally applied to 
the generation of a variety of nonclassical states of many 
photons - in particular multipartite entangled states. 

Here, we propose a new paradigm of entanglement gen- 
eration, which can be implemented in a range of com pac t 
solid-state systems including coupled micropillars 
coupl ed mes as 
ities 
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23| and coupled photonic crystal cav- 
By accurate theoretical modelling 
of the open quantum system, we show that continuous 
variable bipartite entanglement can be generated by an 
array of three weakly nonlinear spatially confined modes, 
linearly coupled via quantum tunnelling. The scheme, il- 
lustrated in Fig. [1^, relies on the quantum interference 
between distinct excitation pathways influenced by the 



sensitivity to small nonlinear shifts of the mode ener- 
gies [iSl ■ This has the advantage of producing degenerate 
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FIG. 1. (color online) a) Illustration of three coupled confined 
optically active modes, with optical pumping of the central 
mode, b) Variation of the entanglement parameter, S\3,, with 
the pump amplitude for different pure dephasing rates, Vp 
(marked on the plot). When the value of Sis is less than unity, 
an entangled state of modes ai and as is present. The dashed 
curve shows the value of 512 evaluated between modes a\ and 
£2. c) Corresponding dependence of the average occupation of 
modes 1 and 3 (solid curves) and of mode 2 (dashed curves), 
d) Variation of Sis with a non-zero detuning between E\ and 
£3- We chose £1+^3 = -0.06meV {Ei = £3 in b and 
c), E2 ~ O.OSmeV. The slight detuning between the modes 
Ei,3 and E2 was found to give the smallest value of 513 for 
fixed J and T by semi-analytic and numerical optimization 
(see supplemental material). Note that J is the largest energy 
scale in the system, such that although the cavities are weakly 
nonlinear they are strongly coupled. 

spatially separated modes, suitable for homodyne detec- 
tion. We demonstrate entanglement by testing the viola- 



tion of inequalities [27 



28| for separable states, and pro- 
vide an intuitive interpretation of how the scheme works. 
Our system is able to reproduce the situation of mode 
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squeezing coupled with beamsplitters [29| in a compact 
microscopic system holding promise for an integrated de- 
vice. The generation mechanism can easily be extended 
to larger arrays of modes, from which multipartite entan- 
glement is expected. 

A general system of three linearly coupled quantum 
boxes is characterised by: energies En{n = 1,2,3); pho- 
ton lifetime fi/F, that we assume equal for the three 
modes; tunneling rate J; and nonlinear energy constant 
U. A near-resonant monochromatic pump drives mode 
2. We assume that the system lies in the weak nonlinear 
regime characterised by [/ < F and U < J 12l4l4|. Note 



able entanglement, continuous variable entanglement is 
characterised by the violation of an inequality (2^, HI] : 



that excellent control over system geometry, energy de- 
tuning and coupling strength has been recently achieved 
experimentally in the case of semiconductor micropil- 
lars [2lj. For near-resonant excitation, higher energy 
modes can be neglected such that each box is described 
by a single mode. Under these assumptions, the system 
is described by the Kerr- Hubbard Hamiltonian [15[: 



n 



(1) 



where a„ are the Bose annihilation operators of the three 
modes and F is the optical pump amplitude. This Hamil- 
tonian is written directly in the rotating frame of the 
pump field, so that is a constant in time and the en- 
ergies En are expressed relative to the optical pump en- 
ergy Hluq. Terms proportional to U describe a Kerr- type 
nonlinearity. These terms are well suited to model the 
Kerr nonlinearity induced by the material (that might 
be enhanced by strong optical confinement) but also a 
resonant nonlinearity due to exciton-exciton interaction 
as e.g. in a confined polariton [l^ system. The quantum 
optical behaviour of our system is fully described using 
the master equation for the density matrix, p: 



..dp 

in — = 
dt 



+ (2"nP^« - flip - pfil) 



(2) 



Two Lindblad type terms account for dissipation at a 
rate F and pure dephasing at a rate Fp, respectively. 
The dissipation is caused by the leakage of photons out 
of the system, while pure dephasing is the result of the 
coupling to a thermal bath [30[ . The latter could be due 
to exciton-phonon scattering in the case of a semiconduc- 
tor structure. Equation [2] can be solved numerically for 
the steady state density matrix using a truncated number 
state basis Q (see supplemental material for details). 

Our aim is to evidence continuous variable entangle- 
ment [2^ between the modes in the first and third quan- 
tum boxes. In analogy to Bell's result for discrete vari- 



1 < S'13 = [pi ~pj,) + V {qi + qs) 



(3) 



where we have defined the amplitude and phase opera- 
tors, pn = {dn + dl) /2 and g„ = (a„ - d^) /(2i), respec- 
tively. The variance of an operator, V{0) = {0'^) — {0)^, 
can be extracted theoretically from the density matrix 
and experimentally measured via homodyne detection. 

For our calculations, we use parameters correspond- 
ing to exciton-polariton boxes although we note that 
the conclusions of our work also apply to several other 
physical implementations. It is well-known how to calcu- 
late the nonlinear interaction strength [ij] and we take 
the value U = 0.012meV in agreement with experimental 
- Ill . A range of coupling strengths are 
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measurements 
possible by varying the separation of the polariton boxes 
and we choose a coupling strength J = 0.5meV, which is 
in agreement with previous theoretical calculations [l9j 
and recent experimental measurements 1211 . The decay 
rate F = 0.044meV was reported in Ref. [34j . 

Figure [T]d shows the dependence of the parameter S13 
on the pump amplitude for a range of values of the pure 
dephasing rate, Fp. For Fp = 0, the black curve shows 
that there is a clear violation of inequality [31 correspond- 
ing to an entanglement of the modes in the first and 
third quantum boxes. In contrast, the modes di and 02 
(or symmetrically 0,2 and 0,3) are not entangled, as ev- 
idenced by the dashed curve showing the value of S12, 
evaluated from Eq. [3] by replacing d^ with 62 • While the 
quantity 5*13 is capable of witnessing entanglement and 
useful given its experimental accessibility, it is important 
to note that it does not fulfill the requirements of a di- 
rect measure of the amount of entanglement (35| . In fact, 
there is no unique, universally accepted, measure of the 
entanglement for our system. 

For increasing dephasing rate, the amount of violation 
decreases and the entanglement is lost at high dephas- 
ing rate. Dephasing rates in semiconductor microcavities 
have been calculated [s^ and measured 37 1 in the range 
of tenths of /xeV. Even for a hypothetical dephasing rate 
an order of magnitude stronger, we still find that the pre- 
dicted violation is sufficient for experimental detection. 

Figure[TJ; shows the corresponding average populations 
of the modes in the signal quantum boxes (solid curves) 
and central box (dashed curves). For small pump am- 
plitudes, corresponding to the linear regime, the popula- 
tions grow according to a power law as expected. Since 
J is large, the largest occupations are those of modes 
Si and 0,3, even though only mode 02 is driven. This 
trend is best understood by expressing the a„ operators 
in terms of eigenmodes of the coupling J. Then, simi- 
larly to the two-mode system [l^ , these eigenmodes are 
driven by the pump in a way that results in destructive 
interference for the occupation of mode 02. Figure [T]i 
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shows the variation of S13 as a function of a finite detun- 
ing between the mode energies Ei and E3. The strong 
resonance at zero detuning is an indication of the un- 
derlying quantum interference mechanism. The level of 
control, required to fabricate a device with such a range 
of detuning to minimize the entanglement parameter S13, 
is achievable in state-of-the-art arrays of semiconductor 
micropillars [2lj |. 

We stress that the reported results are also of sig- 
nificance in several other systems. Since the Jaynes- 
Cummin gs m odel can be linked to an effective Kerr non- 
linearity [38[ , Eq. [1] is also applicable to quantum dots 



embedded in nanocavities and circuit QED systems [39[, 
where the value of U is related to the cooperativity pa- 
rameter. In addition, the value of U has been recently 
evaluated in passive nanocavities [40| . which represent a 
particularly promising system given the low decay and 
dephasing rates. Values of J and Tp suitable for the 
present proposal have also been measured for photonic 
crystal nanostructures. As an example, the coupling of 
nanocavities has been recently studied in Ref. 4l| and 
an upper bound to dephasing rates in quantum dots of 
1/icV has been experimentally established (42| . 

In order to better understand the origin and the na- 
ture of the observed entanglement, we carry out an ap- 
proximate analysis by expanding the quantum state on 
a truncated set of photon number states and solving the 
time-dependent Schrodinger equation for this state. This 
approach does not include the effect of (Lindblad type) 
dissipation and pure dephasing, and is expected to give 
an upper bound to the violation of inequality [H The 
expansion reads: 



iri2n3 

\n1n2n3) 



(4) 



where the basis vectors: 



\n1n2n3) 



000)/^ 



ni\n2\n3\ (5) 



represent states with ni, 71,2 and particles in modes 1, 
2 and 3, respectively. For the analysis, expansion 2] has 
to be truncated to a maximum occupation, N = X)"-*- 
The first ten states, used in expansion 21 arc depicted 
schematically in Fig. [5J together with their couplings 
caused by the pump and tunnelling terms in the Hamil- 
tonian. The states containing two quanta in the same 
mode experience slight energy shifts by an amount 2U 
above the bare energy levels (shown in gray) due to the 
local nonlinear interactions. 

The Schrodinger equation, ihd\%lj) / dt = Hlip), can be 
solved iteratively under the assumption of small occu- 
pations (see the supplemental material for more details) 
for the steady state (including the effect of particle loss). 
The coefficients Cmn^n^ are then calculated and shown 
in Fig. |3] for the cases with ([/ ^ with green/light gray 
bars) and without ([/ = with red/dark gray bars) non- 
linearity. 
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FIG. 2. (color online) Energy level diagram of the first ten 
photon number (Fock) states. The states are labelled by the 
number of photons in each of the three modes. The transitions 
between modes caused by optical pumping and quantum tun- 
nelling are illustrated by the red (dark gray) and green (light 
gray) arrows, respectively. The nonlinear shift of states con- 
taining two particles can be seen by their difference with the 
energy levels calculated in the limit (7 = 0, which are shown 
in grey. 
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FIG. 3. (color online) Quantum state of the system expanded 
onto the particle number states. The values of |C„in2"3p 
are plotted in the linear (red/dark gray bars) and nonlinear 
(green/light gray bars) regimes. Note that not all the ba- 
sis states are shown, since due to symmetry |Cnin2Ti3l — 
|Cn3Ji2"iP- The states with 712 = have the highest occu- 
pations and their relative occupations are exactly given by 
the binomial coefficients in the linear case (dashed lines and 
black points). The parameters were the same as in Fig. [TId 
with F = O.SmeV and F = 0.044meV. 



In accordance with Fig. [T];, we observe that the quan- 
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turn state is in general characterised by very low occu- 
pancy of mode 2. Each photon that is initially injected 
in this mode, tunnels to modes 1 and 3. This behaviour 
can be easily understood in the linear case {U = 0), for 
which the Hamiltonian can be diagonalised exactly. In 
this case, the Schrodinger equation shows that only the 
mode generated by the operator (oi + 03)^ is effectively 
driven by the pump, thus giving rise to a fully separa- 
ble quantum state, expressed as a linear combination of 
states ((Zi -|- 03)^''^ |000) at varying occupancy N. Conse- 
quently, the relative weights of the coefficients Cnin2n3, 
for each given value of the total occupancy N, are ex- 
actly given by binomial coefficients, as shown by the 
dashed lines in Fig. [3] In the nonlinear regime, the sys- 
tem changes to a state characterised by the green (light 
gray) bars in Fig. [31 where it is clear that states con- 
taining particles in both modes 1 and 3 (e.g., |101)) are 
suppressed with respect to the linear case, while those 
with all particles in the same mode are enhanced (e.g., 
1 200)). This result is a consequence of the nonlinear shift 
when photons occupy the same box, which has an ef- 
fect on the quantum interference of possible time evo- 
lution paths in the Fock basis. As an example, within 
the manifold of states with N = 2 occupancy, the cou- 
plings of the states |110) and |011) to the states |200) 
and 1 002), respectively, change the phase of any time- 
evolution path passing through those states. If we con- 
sider the system initially in the state |011) for example, 
then to reach the state |101) two possible options are 
clear from Fig. [2j the direct path |011) i-> |101) or the 
path |011) ^ |020) ^ |110) ^ |101). The relative quan- 
tum phase of the two paths is affected by coupling to the 
state 1 200), which in presence of the small nonlinear shift 
of this state induces a destructive interference of the two 
paths and suppression of the state |101). We arc left with 
a situation where the detection of at least one photon in 
either signal mode, 1 or 3, grants that no photon will be 
detected in the other mode. This result solely depends on 
the nonlinearity in modes 1 and 3. We have verified that 
the parameter 6*13 experiences negligible change when the 
nonlinearity in mode 2 is removed. 

In summary, arrays of coupled photonic modes are able 
to display striking quantum correlations despite their 
modest nonlinearity in the low occupation limit. This 
allows continuous variable entanglement to be generated 
between degenerate spatially separated modes that are 
coupled via quantum tunnelling, in a way that is robust 
to typical decoherence rates in these systems. The set of 
three coupled modes here described serves as a building 
block that can be repeated on an array of modes with 
appropriate topology, which could be further controlled 
using electric or magnetic fields [i^. This sets a viable 
paradigm for the generation of multiparty entanglement 
in arrays of quantum boxes on a single device. 
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